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0^ ■ ABSTRACT 

(^ 

C i The distance to the Large Magellanic Cloud (LMC) can be directly 

determined by measuring three of its properties, its radial- velocity field, 
its mean proper motion, and the position angle 0ph of its photometric 
line of nodes. Statistical errors of ~ 2% are feasible based on proper 
motions obtained with any of several proposed astrometry satellites, 
CS ! the first possibility being the Full-Sky Astrometric Mapping Explorer 

(FAME). The largest source of systematic error is likely to be in the 
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O ! determination of 0ph. I suggest two independent methods to measure 

Q^ ! 0ph, one based on counts of clump giants and the other on photometry 

of clump giants. 1 briefly discuss a variety of methods to test for other 

^i sources of systematic errors. 
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1. Introduction 

The distance to tlie Large Magellanic Cloud (LMC) plays a crucial role in 
the extragalactic distance scale. The relation between log-period and apparent 
magnitude of LMC Cepheids is quite well determined. If an LMC distance c^lmc and 
mean LMC-Cepheid reddening are assumed, then the Cepheid period-luminosity 
relation is effectively calibrated. The distance to external galaxies harboring 
Cepheids can then be determined by comparing their observed fluxes to those of 
LMC Cepheids at the same period, and by taking account of the differences in 
reddening which are determined from the differences in color between the target 
Cepheids and those in the LMC. A variety of secondary distance indicators have 
been calibrated in this fashion. While a decade ago, the Hubble constant Hq derived 
from these measurements ranged over a factor of two depending strongly on both 
the author and the method, a major observing campaign with the Hubble Space 
Telescope (HST) has dramatically narrowed this conflict. For example, Saha et al. 
(1999) recently find i^o^LMc/SO kpc = 60 ± 2 (internal) kms^^ Mpc^^ while Madore 
et al. 1999 find ifo'^LMc/50 kpc = 72 ± 3 (random) ± 5 (systematic) kms""*^ Mpc^-*^. 

By constrast, the disagreements over c^lmc have not narrowed at all over the past 
decade. The primary methods for measuring c^lmc use "standard candles", objects 
whose luminosity is presumed to be fixed or to depend only on distance-independent 
observables such as period, metallicity, etc. Their absolute magnitudes must be 
calibrated locally. Two major standard candles that have been used to measure c^lmc 
are Cepheids and RR Lyrae stars. Three recent determinations, all from Hipparcos- 
based calibrations of these standard candles, illustrate the range of o^lmc estimates: 
Feast & Catchpole (1997) find c^lmc = 55.0 ± 2.5 kpc based on trigonometric 
parallaxes of Cepheids; Gould & Popowski (1998) find c/lmc = 45.1 ± 2.7kpc based 
on statistical parallax of RR Lyraes [and assuming Vo{RR,LMC) = 18.98, Walker 
1992]; and Gratton et al. (1997) finds c/lmc = 52.1 ± 1.7kpc from an RR Lyrae 
calibration based on fitting globular cluster main sequences to Hipparcos subdwarfs, 
while Reid (1997) finds a slightly longer distance based on the same technique. At 
present it is not known if these discrepancies are due to undetected systematic errors 
in the various techniques or to non-standardness in one or more of the "standard 
candles", or both. It is unlikely that that the differences are merely statistical 
fiuctuations. For example, Popowski & Gould (1999) review a variety of methods for 
calibrating RR Lyrae stars whose results disagree by substantially more than their 
statistical errors. 



Of course, one would prefer to eliminate the distance ladder altogether and 
simply obtain a direct measurement of ^lmc There are two possible paths to a direct 
distance measurement: trigonometric parallax and kinematic methods. The parallax 
of the LMC is TT ~ 20 /xas. The Space Interferometry Mission (SIM) should be able to 
make individual astrometric measurements accurate to ~ 8;uas, and could perhaps 
achieve a^ ~ 2 /las given a sufficient number of observations. This limit is set by the 
precision of the SIM "grid-star" solution. Hence it cannot be significantly improved 
upon by making measurements of several LMC stars, since these lie in the same 
field. While such a ~ 10% meausurement would certainly be of interest, it would not 
by itself clearly distinguish among the various competing distance estimates. 

Panageia et al. (1991) made the first kinematic measurement of c^lmc by 
comparing the light travel time accross the ring around Supernova 1987A with 
its angular diameter as measured by HST. They found /ign = 18.55 ± 0.13 where 
/isn is the distance modulus of Supernova 1987A. Gould (1995b) reanalyzed these 
data and obtained /isn < 18.350 ± 0.035. Sonneborne et al. (1996) rereduced the 
original light-curve data and found /ign = 18.43 ± 0.10. Gould & Uza (1998) then 
reanalyzed these rereduced data and obtained /isn < 18.372 ± 0.035 if the ring 
were assumed circular, but /isn < 18.44 ± 0.05 if it were assumed elliptical (as 
some evidence suggests). Finally, Panagia (1998), using the same data, but arguing 
that the effective radius of the ring had grown between the time of the light echo 
measurements and those of the angular size of the ring, found /isn = 18.55 ± 0.05. 
In brief, there remains controversy over the interpretion of the data at the ~ 10% 
level in distance. Since the event itself was unique and the measurements will 
never be repeated, it seems unlikely that this conflict will be resolved to everyone's 
satisfaction. 

Here I propose to use the radial-velocity gradient method to measure the 
distance to the LMC. This method has been used in the past to measure the 
distance to the Hyades (Detweiler et al. 1984; Gunn et al. 1988) and the Pleiades 
(Narayanan & Gould 1999). When applied to the LMC, the method has some unique 
characteristics relative to previous applications. This is in part because the LMC is 
a cold system supported primarily by rotation while the Hyades and Pleiades are 
supported by pressure, and in part because the LMC is two orders of magnitude 
farther away. 

Consider a cold disk rotating at a projected angular rate Q{R), and moving with 
a systemic proper motion fi (and so transverse velocity V_|_ = nd). li fi <^ Q, then 



the locus of extrema in the radial velocity field will coincide with the photometric 
major axis. That is, the kinematic and photometric lines of nodes will be aligned. 
However, the transverse velocity V_l of any system induces a gradient in the radial 
velocities because the radial vector that is dotted into the velocity to form the radial 
velocity changes direction across the system. Thus the observed gradient Vvr will 
be displaced from the that due to internal rotation alone by V^. If the direction of 
the photometric line of nodes is known, and if fi is measured (so that the direction 
of V_L is also known), then it is straight forward to solve for the magnitude of 
V_|_. The distance is then simply d = Vj_/fi. The LMC is sufficiently close and is 
moving sufficiently rapidly, that the kinematic line of nodes is displaced from the 
photometric line of nodes by about 25°. The interpretation of this displacement 
could be clouded by uncertainty about how well the LMC conforms to the ideal of a 
flat axisymmetric system that I use to model the data. After I present the method 
and derive the statistical uncertainties, I briefly discuss how the measurement could 
be corrupted by systematic deviations from this ideal, and I indicate some methods 
to check for such systematic effects. 

To illustrate the method, I will assume the use of astrometry data such as 
would be obtained by the Full-Sky Astrometric Explorer (FAME), a proposed Midex 
mission. As I will discuss, the method could also be applied to data from SIM or the 
Global Astrometric Interferometer for Astrophysics (GAIA). 



2. The Method 

Consider a stellar system whose physical size is small compared to its distance d. 
Let the space velocity of the system be V and let the space motion of an individual 
star in the sysetm be Vj. I then write 



V + u,- + (5v, 



where Uj is the mean internal systemic motion of the stars in the system (due, e.g., 
to rotation) at the projected position of star i, and (5vj is the peculiar motion of star 
i relative to this systemic motion. The radial-velocities are therefore given by, 

Vr,i = rij • V + Ur,i + 6Vr,i (2) 

where iij is the unit vector in the direction of star i, Ur^i = n, ■ Uj, and 6vr,i = rij ■ (5vj. 
I assume that the radial-velocity residuals 6vr,i are randomly distributed with 



dispersion a^- I also assume that the internal systemic motion Uj is known. In fact, 
determining u is not trivial, but I ignore this problem here and return to it in §§ 3.2, 
3.3, and 4. Then the radial-velocity gradient with respect to angular position on the 
sky is given by 

VVr = Vx + VUr (3) 

where V^ = V — no (no • V) is the transverse velocity of the center of the system, 
and no is the direction vector pointing to this center. Since Vfr is a vector, the 
errors are properly described by a covariance matrix, Cxy This is given by 



c = b ^ bki = -^ 
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where = {6x, Oy) is the angular position of star i relative to the center, and A^^ is 
the number of radial-velocity measurements. However, for simplicity, I will consider 
stars distributed uniformly over a circular area of radius A9. In this case, the error 
in each component of Vvr (or equivalently V^, since Vwr is assumed known) is 

See also Narayanan & Gould (1999). 

Suppose now that the proper motion /i of the system is measured with error a^ 
The distance and distance error are then. 



jji \ d J \V±J V/i 

That is, the fractional distance error is limited by the larger of the errors in the 
transverse velocity and the proper motion. 



3. Application to the LMC 
3.1. Naive 

Substituting values appropriate for the LMC into equation (|[), I obtain 

V< °UOV 24kms-iV325kms-V V 4° / ^^ 



where I have chosen a dispersion characteristic of carbon stars (Cowley & 
Hartwick 1991) and the estimate of the transverse velocity from the proper-motion 
measurement of Jones, Klemola, & Lin (1994). Hence, good statistical precision 
is possible provided that a large sample of stars is available. Note that the 
measurement errors are not important provided that they are well below dispersion. 
Since a.,; ~ 5kms~^ errors are not difficult to achieve for LMC carbon stars, it is 
feasible to obtain a very large sample such as is envisaged in equation (|^). 

While the proper motion of the LMC is only crudely known today (Jones et al. 
1994; Kroupa & Bastian 1997), it could be measured to very high precision with any 
of a number of proposed astrometry satellites including FAME, SIM, and GAIA. 
For definiteness, I will focus on the capabilities of FAME which has the earliest 
possibility of launch. I find from the USNO-A2.0 catalog (Monet 1998), that there 
are a total of 21,900 stars with 13 < l^ < 15 within A6 = 4° of the center of the 
LMC at {l,b) = (280.5, -32.9), where I estimate V = {B + R)/2. Of these, about 
13,300 are foreground Galactic stars as judged from counts in three similar circles at 
{l,b) = (280.5, +32.9), and {l,b) = (79.5, ±32.9). This leaves N^ ~ 8,600 stars in 
the LMC. The dispersions of LMC stars in the transverse directions are unknown, 
but based on what is known of disk kinematics in the Galaxy, it is plausible to 
assume that they are ~ 50% higher than the vertical dispersion, or a± ~ 35kms^^. 
Hence, if the proper motions of these stars could be measured to better than 
o'±/dhMC ~ 150/iasyr~^, and if the internal systemic motions u are again assumed 
known (see § 4), then the precision of the LMC proper motion would be given by 



Arr^/";± = 0.l% 



/i V^ ' V8600y 35kms-iV325kms-i 

where N^ is the number of proper- motion measurements. In fact, FAME probably 
cannot achieve quite this precision at V = 15, but should come within a factor of 
2 (Horner et al. 1998) and so easily achieve a^j^ln <^ 1% or a^ <^ 10/iasyr^^ The 
present rotational precision of the extra-galactic reference frame is a^ ~ 5/iasyr~^. 
However, the FAME astrometric frame will probably be accurate only to within 
cr^ ~ 25/iasyr~^ (assuming 100 QSOs with V ^ 15 and hence with mean proper 
motion errors of 250/xasyr~^). The FAME frame will be fixed substantially better 
by SIM. In brief, the proper-motion measurement error can probably be reduced to 
about 2% with FAME alone and substantially less by combining FAME and SIM. 



3.2. Degeneracy 

However, the potentially fatal flaw in this method is that u is not known (as has 
been assumed so far) and must be determined from the same kinematic data that 
are used to derive the distance measurement. As is well known from the classical 
application of the radial-velocity gradient method to the Hyades (Detweiler et al. 
1984; Gunn et al. 1988) and the Pleiades (Narayanan & Gould 1999), if the cluster 
were undergoing solid-body rotation u = H x r, this would produce a radial-velocity 
gradient 

Vur = (no X Q,)d. (9) 

Here r is the 3-space position of a star relative to the cluster center. This gradient 
is indistinguishable from the gradient produced by a transverse velocity and so, if 
unrecognized, would corrupt the distance measurement given by equation (^. In 
the case of clusters, one normally simply assumes that the cluster is not rotating. 
However, one can check this assumption by comparing the directions of the radial 
velocity gradient and the proper motion. If these differ, the cause might be rotation 
(or systematic errors). If they are the same, then either the cluster is not rotating, or 
its rotation happens to be perfectly aligned with its proper motion (within statistical 
errors). 

The situation is similar for the LMC but is somewhat more complicated because 
the LMC is rotating. While the rotation is not solid body, it can be reasonably 
approximated as such in the inner 2°.5. To the extent the rotation is solid-body, one 
measures a gradient 

Vfr = V_L + fixC^LMC, fix = Hq X 11, (10) 

and from this measurement alone, has no idea how to separate the two components. 
If, for example, one ignored the transverse motion, one would interpret the gradient 
as due entirely to rotation and would therefore misjudge the amplitude of rotation. 
One would misjudge its orientation as well to the extent that V^ does not happen 
to lie parallel to Hx- 



3.3. Breaking the Degeneracy 

However, for a disk rotating about its axis of symmetry, Q^ should be aligned 
with the apparent major axis of the system, i.e., the photometric line of nodes. This 



provides some information with which to break the degeneracy. 

These effects were first investigated when Feitzinger et al. (1977) reanalyzed 
earher kinematic data. They noted that the kinematic hne of nodes (locus of extrema 
in radial velocity) was displaced by ~ 20° from the photometric line of nodes (major 
axis of the surface-brightness profile) at position angle 0ph = —10°. They assumed 
that this displacement was caused by transverse motion in the direction 0^ = 110° 
(i.e., the direction of the Magellanic stream) and then solved for the amplitude of 
this motion V± ~ 275kms~^. Subsequently, several other workers applied a similar 
procedure to various stellar samples and obtained various results (Rohlfs et al. 1984; 
Meatheringham et al. 1988; Hughes, Wood, & Reid 1991). Note that this approach 
to breaking the degeneracy requires two pieces of information in addition to the 
kinematic data: first the position angle of the photometric line of nodes 0ph, and 
second the direction of LMC motion (p^. 

However, if the proper motion is measured (which is necessary in any case to 
determine the distance through eq. H), one already knows 0^. From equation ([10|) , 
the three vectors Vwr, V_|_, and HxC^lmc form a triangle, so by the law of sines, 

VI = ""[^-M iv-^,!. (11) 

sm{(pf, - 0ph) 

where 0v is the observed position angle of the kinematic line of nodes. The 
quantities on the right-hand side of equation (|T^) are all observables. Assuming that 
the errors in the measurements of fi and Vvr are isotropic, so that o"(0^) = cr^/fi 
and cr(0v) = crv/|Vfr|, one can evaluate the error in ^lmc = ^±/yU by taking the 
derivatives of equation ( p^T]) with respect to the various parameters. I find. 
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^ph, (12) 



_sm(^(/)v - (pph) sm[(pf, - (pph) 
where Uph is the error in the determination of </)ph. 

Equation (|12D differs from its naive relative, equation (|^), in two ways. First, the 



entire error in equation (^ is now multiplied by a factor csc(0^ — 0ph)- Second, there 
is a new term which is related to the uncertainty in the photometric position angle. 
To understand the importance of these changes, I first introduce representative 
values of the parameters. I choose 0^ = 97° from the proper motion meausrement 
of Jones et al. (1994), fix = 12kms~^ kpc~^, and Vj_ = 325kms~^. Together, these 
imply 0v = 14°, thus (0^ - 0ph) = 107°, (0v - 0ph) = 24°, and (0^ - 0v) = 83°. 



The fact that cj)^ and (/)ph are almost at right angles implies that the csc^(0^ — 0ph) 
term in equation (p!2D is essentially unity. However, since fixdhuc ^ Kl, the 
radial-velocity gradient due to V^ is a relatively minor perturbation on the gradient 
due to internal motion, and so 0v is not much different from 0ph. Hence, the 
factor sin(0v ~ 0ph) = 0.41 in the denominator of the last term is relatively small. 
This means that 0ph must be measured quite accurately if one wants a precise 
measurement of the effect of the transverse velocity. Specifically, the last term in 
equation (|I^) is (2.55o"ph)^. At distances from the center ^ 2°.5, the rotation curve 
tends to flatten, and so V_|_ becomes a larger relative perturbation causing (0v ~ 0ph) 
to grow and thus making the measurement somewhat easier. Nevertheless, imprecise 
knowledge of 0ph is likely to be a major limitation of the method. 



4. Measurement of 

To achieve 2% precision in (Jd/d (which generally seems feasible from the 
standpoint of the Vf r and jj, measurements) would require measuring the position 
angle to o"ph ~ 0°.4, or 0.008 radians. It is difficult to believe that this can be 
achieved using surface photometry alone. Recall, that one is not actually interested 
in the best fit to the major axis of the isophotes. Rather, one wants to know the 
position angle of the line that crosses the plane of the sky. Certainly star formation, 
dust, etc corrupt the surface-brightness profile too much to extract information at 
this level of precision. It should be possible to make a more accurate assessment 
of 0ph using star counts particularly of clump giants. Using the method of Gould 
(1995a) one may show that this technique can determine 0ph with precision 



O'ph 
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where N^g is the number of clump giants, i is the inclination of the disk, F{R) is 
the (assumed axially symmetric) radial profile of the LMC disk, and where I have 
assumed i = 30°, and {{dlnF/dlnRY) = 6, which is valid for an exponential disk. 

However, clump giants provide another, independent route to the measurement 
of the position angle. Clump giants have a dispersion in / band of only Ccg = 0.15 
mag (Udalski et al. 1998). The stars on the near side should therefore be brighter 
than those on the far side by a significant fraction of this dispersion. Averaging this 
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effect over tlie wliole disk, I find tliat (jph can be determined to a precision 



wfiere I liave assumed an exponential scale length of a = 1°.7 (Feitzinger et al. 
1977), so that {R'^)/d1^Q = 6a^. The challenges to actually carrying out such a 
measurement would be formidable. Just maintaining a constant photometric zero 
point at the level of o"cg/^cg^ ~ ^0"^ '^^S ^^^^ fields separated by ~ 10° would be 
difficult. In addition, one would have to correct for differential reddening, probably 
from the clump giant colors, but to do so would require an accurate estimate of 
E{V — I)/Aj. This could be made empirically by looking at the correlation between 
V — I and / at fixed position but might not be easy. 

In principle, it is also possible to measure 0ph from the transverse velocity 
field measured from the proper motions. In practice, however, the errors in this 
determination are too large for it to be useful. Note that the internal transverse 
motions do not increase the uncertainty in /x. The uncertainty in the transverse 
velocity field at any particular point is much smaller than either the dispersion or 
the measurement error, and there is no uncertainty in the mean internal motion 
averaged over the whole population: the mean internal motion is zero. 



5. Discussion 

1 have outlined how the radial-velocity gradient method could be applied to 
measure c^lmc with statistical errors of 2% or less. Of course, as in most distance 
measurements beyond the solar neighborhood, the largest potential source of errors 
is systematics. Examples of effects that would generate such systematic errors are 
non-circular motions and/or warps in the LMC disk and contamination by material 
along the line of sight. For example, Weinberg (1999) has recently shown that 
resonant interactions between the Milky Way and the LMC can profoundly disturb 
the LMC disk. 

However, given the mass of data required to make the measurements, it should 
be possible to conduct many tests for systematics. For example, non-circular motions 
would affect both the radial- velocity gradient and the orientation of the photometric 
line of nodes. The latter would have a larger impact on the distance simply because 
its coefficient in equation (|T^) is ~ 2.4 times larger. Such motions should be revealed 
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in the comparison of the clump-giant star-count and photometric methods for 
measuring 0ph: the star-count method would be affected by non-circular motions 
while the photometric method would not. Both warps and non-circular motions 
could be tested by comparing the radial-velocity field with the transverse-velocity 
field obtained from proper motions. Similarly, it is possible to search the radial 
velocities for evidence of unassociated material along the LMC line of sight (Graff et 
al. 1999). 

While an all-out serach for systematic effects probably requires the full data 
set, substantial initial investitgations can be made with existing phtometric catalogs 
or with radial velocity studies now underway (e.g., Suntzeff, Schommer, & Hardy 
1999). 

I have estimated that FAME will obtain 8600 proper motions with a mean 
precision of 250/xasyr~^. If FAME is not launched, what are the prospects for 
matching this performance? Clearly GAIA, which is also a survey mission but with 
much higher precision and fainter magnitude limits could easily meet this standard. 
However, given its later launch data and larger analysis time, GAIA would require 
an additional decade to produce results. SIM certainly has the capability to make 
these measurements, but whether it would make so extensive a survey is open to 
question. Recall that the proper-motion measurements need only be a factor of a 
few better than the internal dispersion (~ 150/iasyr~^). For y ~ 15 stars, SIM 
could do an order of magnitude better than this in 1 minute. Allowing another 
minute for pointing and assuming a total of 4 position measurements per star, 8600 
proper-motion measurements would require about 1000 hours. From equation (|^) 
only ~ 100 stars would be needed to measure c/lmc to ~ 1%, and this could be done 
in only about 10 hours. In this case, however, one would lose much of the ability to 
check for systematics from a comparison of the radial-velocity and proper-motion 
fields. In brief, FAME is the instrument of choice to make the proper-motion 
measurement. 
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